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Abstract 

We address the dynamics of damped collective modes in terms of first and second mo- 
ments. The modes are introduced in a self-consistent fashion with the help of a suitable 
application of linear response theory. Quantum effects in the fluctuations are governed 
by diffusion coefficients D^ u . The latter are obtained through a fluctuation dissipation 
theorem generalized to allow for a treatment of unstable modes. Numerical evaluations 
of the D^ u are presented. We discuss briefly how this picture may be used to describe 
global motion within a locally harmonic approximation. Relations to other methods are 
discussed, like "dissipative tunneling", RPA at finite temperature and generalizations 
of the " Static Path Approximation" . 
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1 Introduction 



The interplay of dissipation and quantum effects has been a long standing challenge to trans- 
port theory, in particular the question of the influence of damping on quantum tunneling. 
Convincing answers have only be found in the 80'ies by clarifying demanding problems of 
principal nature with the help of functional integrals, see e.g. 0-0. In these approaches 
one first derives effective actions for global motion from which one may then deduce expres- 
sions for the decay rate of a meta-stable state. In contrast to this procedure, we like to look at 
the problem from the point of view of transport equations such as the one of Kramers || but 
extended to include quantum effects. By its very nature, such a transport equation reflects 
local properties of collective motion as displayed through the propagation in collective phase 
space. In the following we like to concentrate on situations where damping cannot simply be 
treated perturbatively. For the case of stable modes, suggestions of possible quantal equations 
have been made for instance in 0j and M. The possibility of generalizing to unstable modes 
has first been described in [[J. 

Different to those papers mentioned previously which deal with dissipative tunneling and 
which exploit functional integrals, we like to base our analysis on systems where collective 
variables are introduced self- consistently. For Hamiltonian dynamics, this implies that in the 
Hamiltonian 

7~L — H intr + H coup i + H co \\ (1) 

for the total system both the "collective" part H co n as well as the coupling term are bound to 
be functionals of the dynamical variables of the intrinsic system. One prime example is the 
electron gas for which D. Bohm and D. Pines have been able to deduce such a Hamiltonian 
in | 10| . Another example is nuclear collective motion for which an adapted version of the 
Bohm-Pines method can be found in | H| , see also [12|]. Although later on we are going to 
concentrate on systems for which such a Hamiltonian (fl|) exists it may be mentioned that our 
method of accounting for self-consistency in the fluctuating force can be generalized to other 
cases, the most prominent example being a situation in which the dynamics is described by 
the Landau- Vlasov equation with (collision term) (see \T3\). For stable modes this method 
may in some sense be considered a generalization of the suggestion made by Landau (see 

EI)- 

Evidently problems like the motion from a meta-stable minimum across a potential barrier 
are manifestly non-linear in at least one of the collective degrees of freedom. Moreover, there 
are systems for which the intrinsic sub-system and its coupling to the collective part change 
along the dynamical path, thus prohibiting to linearize this coupling in the bath variables. 
In such a situation, as it is given for example in nuclear fission, the technique of functional 
integration is not useful anymore as the internal degrees of freedom can no longer be integrated 
out. It has been argued that it might be possible to handle such systems by a locally harmonic 



approximation (LHA ) (see e.g. [16|, |fL7|l , [|12|| ) . In this method it is not necessary to restrict the 
dynamics of the internal degrees of freedom to a system of (perturbed) oscillators. Certainly, 
as an essential intermediate step, there will appear the dynamics in linearized version for 
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the collective degrees of freedom. As already demonstrated by Kramers, in lucky cases it 
may even suffice to have information on the dynamics at the potential minimum and at the 
barrier, the latter mode being unstable. In this paper we shall address specifically these 
problems and study the influence of dissipation on average motion as well as on the dynamics 
of the second moments (fluctuations). This will be achieved within the framework of linear 
response theory. The latter does not only allow one to treat average motion (in a linearized 
fashion), it opens an ideal playground for deducing the diffusion coefficients which govern 
the dynamics of the second moments. To be able to treat instabilities as well, it will turn 
out necessary to perform suitable analytical continuations. Manipulations of the kind which 



will be needed here have been used before in |Iq| , but, as mentioned, not in connection 
with a transport equation and, hence, diffusion coefficients. As we shall see the numerical 
computations of these coefficients will very naturally hint to limits of the LHA at smaller 
temperatures. 



2 Brief review of average motion 

In this paper we would like ourselves to restrict to a situation in which the dynamics of interest 
can be described by the one quantity F. For average motion this means looking at the time 
evolution of the expectation value (F) t . To start with a simple situation let us suppose we 
are given a Hamiltonian of the type 

H {2) = H + k -FF (2) 

with some coupling constant k, without bothering for the moment how it has come about. 
Later on we are going to describe how this form may derived from more general considerations 
and how the k may be deduced from static properties of the system, in which case it will 
depend on the actual state of the intrinsic system. Let us suppose further that we are 
interested in harmonic vibrations. As may be well known, their frequencies and strengths can 
be deduced from the following response function 

obeying the following definition 

6{F) U = -XcollH/extM (4) 

The S(F) UJ measures the Fourier transform of the deviation of (F) t from the equilibrium value 
(F)o, as it follows from the response of the system in linear order to an external field f ext (uj), 
with the coupling being given by 5H = f ex _tF. The x{ u ) is the Fourier transform of the 
response function defined by 

X(t - s) = &(t - s) Z -tT (p2 s [F 7 (t),F'( S )]) (5) 
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with the /)q S representing the equilibrium density operator associated to the Hamiltonian H . 
Notice please that this function is governed entirely by intrinsic properties, different to the 
Xcoii(a>) which contains information on collective motion. 

We should expect average dynamics to be closely related to the mean field approximation. 
Indeed, for the latter one readily verifies the following Hamiltonian 

H$ = H + (Q- Qo)F (6) 

where the abbreviation 

k(F) t = Q-Q (7) 

has been used. One may now simply apply the Clausius-Mosotti construction to derive the 
form (^) from @ and ([/]). 

For the following it will show convenient to add to the Hamiltonian H^J of (Q) a c-number 
term and thus to introduce a H(xi,pi, Q) like 

H(xi,pi, Q) = Qo) + (Q~ Q )F(x l ,p u Qo) + \{Q~ Qo) 2 ( d2H{ QQ2 l,Q) )Q ■ (8) 

Identifying H(xi,pi, Q ) = H and the expectation value on the very right with 

_, / d 2 H(x i ,p l ,Q) ^ 



one sees that the total energy of the system E tot = (H^J) — -^(F)(F) can simply be expressed 
as the expectation value (H(xi,pi,Q)) of the one-body Hamiltonian H(x{,Pi,Q). 

This discussion has demonstrated the intimate relation between a Hamiltonian of the form 
(f|) and the one of (|2|) involving a schematic two-body interaction. On the level of the mean 
field approximation the latter becomes equivalent to the former. However, the of (0) 
allows one to treat fluctuations in F, or what will turn out more convenient in Q — Q , about 
its average value. Indeed, later on we will be exploiting exactly this feature when we are 
going to quantize the "collective variable" Q. Conversely, as often mean field approximations 
are more easily accessible numerically than solutions of the full Schrodinger or Heisenberg 
equations (of the many body problem), one may start from such a H(xi,Pi, Q) and construct 
the effective two body interaction appearing in H^ 2 \ 

A prominent example from nuclear physics is given when one tries to simulate the mean 
field approximation via the deformed shell model. In that case the Q is introduced as a 
c-number variable specifying the shape of the nucleus. The latter then appears both in the 
single particle potential as well as in the liquid drop energy. To account for the latter is 
necessary because otherwise one would not be able to calculate reliably the total energy of 
the system. This can be done with the help of the Strutinsky renormalization. How one 
may then proceed to construct a Hamiltonian of the type H(xi,Pi,Q) has been worked out 
in detail in |2(| . 



It is more than conceivable that one might start from even more microscopic theories like 
the one of Hartree-Fock involving effective forces or the Relativistic Mean Field theory. An 
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equation like ([/]) would then appear as a constraint for possible variations in the deformation 
of the self-consistent field. After appropriate linearizations and after applying time dependent 
perturbation theory in suitable manor one would end up in a form like @, provided that one 
is interested in the change of the expectation value of the same operator F which specifies 
the variation in the density. In general the response "function" would have tensor character, 
of course. 

In this context one may mention examples worked out for Landau Fermi liquid theory 
21] , fl3| . Looking at general modes of the single particle density in phase space for the 



in 



associated response functions Xcoii(^) forms similar to the one given in (|3|) have been deduced. 
There it might only happen that the " coupling constant" k may depend on the wave vectors 
of the various modes. 



2.1 Extension to finite excitations and large scale motion 

In the following we would like to concentrate on the example for which a Hamiltonian 
H(xi,pi,Q) is given, which depends parametrically on a c-number Q(t) representing the col- 
lective degree of freedom. It is this example which is developed farthest theoretically and 
which still can be handled numerically even for such complex situations as nuclear fission 
[ ^| , As described in various papers (see e.g. |^4| and 0], in particular), with such a 



Hamiltonian it is possible to derive an intrinsically consistent theory for collective motion of 
large scale provided there is a clear separation of time scales. 

Let r co ii be a relevant measure for a collective time scale and r the one for the remaining 
degrees of freedom, henceforth called the intrinsic ones. For time lapses St with t < St 
t co \\ one may describe collective motion within a (locally) harmonic approximation. This 
is achieved by exploiting the expansion (||) and identifying simultaneously the unperturbed 
density operator p° s as the one corresponding to the Hamiltonian H(xi,pi,Qo). This density 
operator may then either be a function of temperature or entropy depending which distribution 
we want to use for parameterizing the equilibrium. For the sake of simplicity, in this paper 
we will choose the canonical one later on, but we may refer to [12| for a discussion of the more 
general case. Thus we may write for the density p^ s = p qs (Qo, To), with the Qo and the T to 
represent those values which the "macroscopic" variables Q,T have during the time interval 
St. 

As shown in [23j] and all the steps discussed above for the genuinely harmonic case 
may then be taken over. We must expect, of course, that all quantities will depend on the 
actual quasi-static state of the system which is being specified by the pair of variables Qo, T . 
This is true in particular for the coupling constant k which can be written in the form 



k 



, _ d 2 E(Q,S 



dQi 



+ x (0) = C(0)+x(0). (10) 



Here E(Q, So) is meant to represent the internal energy of the system, with the entropy Sq 
being understood to be calculated for given Q ,T . (It is possible, of course, to rewrite this 
formula in terms of derivatives of the free energy). The x(0) stands for the static response and 
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the C(0) is introduced as a short hand notation for the local static stiffness. To derive this 
relation, together with the form (|3]) for the collective response function, one needs to apply 
similar arguments as mentioned before but has to care for possible changes in the quasi-static 
properties of the system. In doing so one may be guided by the fact that within the harmonic 
approximation (for the motion in Q) entropy can be assumed to be constant (to leading 
order), as it will change with velocity only quadratically. Using this hypothesis, in formal 
sense the derivation of the form of the collective response Xcoii(^) follows identically the one 
given above and which ended in (§), provided the adiabatic susceptibility % ad is identical to 
the static response, x ad = f° r a general discussion of the relevance of this property see 
2IJ and (T| . 



The existence of equilibrium for the nucleonic degrees of freedom requires the presence of 
residual interactions. However, it would be beyond the scope of any realistic theory to allow 
for both a detailed treatment of such interactions as well as numerical applications, say for a 
calculation of the transport coefficients which such a theory aims at. It has been suggested, 
therefore, to acount for that by dressing the particles and holes (of the deformed shell model) 
by a self energy whose imaginary part is given by 

r(u>,T) = — - \ w 11 
r 1 + ^ [{huj - /i) 2 + tHt 2 j 

with /i being the chemical potential. It is this imaginary part which finally will be responsible 
for the presence of macroscopic damping. It would out of scope of the present paper to explain 



more details; instead, the reader is referred to | 12fl . In numerical applications the following 
values have been used for the two parameters: Tq = 33MeV and c = 20MeV. 

With respect to the choice of Qq, we should like to note that for the general situation 
away from the potential minimum, the Q in (^) is to be replaced by the center Q m of the 
oscillator approximating the "true" static energy E(Q,S ) in the neighborhood of the Q 



of (Rp (see e.g.JT2j]). Furthermore, we want to mention that the discussion to come shall be 



restricted to situations of negative coupling constants, which for the case of an instability 
will imply the | C(0) | to be smaller than x(0); see below. Looking back to the two body 
interaction introduced in (]|) this is to say that we like to restrict ourselves to effective forces 
being attractive. For the nuclear case this is known to correspond to motion with neutrons 
and protons in phase. 

2.2 Reduction to single modes 

The poles of Xcoii(^) define possible excitations of the system, more precisely those which 
involve dynamics in the quantity (F) and which are caused by that perturbation SH which 
is proportional to the F itself. Considering all of them would in the end imply strict non- 
Mar kovian behavior in (F) t . Eventually the strength distribution is dominated by only one 
mode, in which case the Xcoii(^) m &y be replaced by the response function of the damped 
oscillator 

XoscM = t\iTF 2 i • r 2 ( 12 ) 

M r to 1 + hoy — w A 
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with M F T = 7 F and M F zu 2 = C F , where 7 F , M F and C F stand for the coefficients of friction, 
inertia and stiffness for the (local) motion in 5(F) t , a feature which is easily verified from (fj). 
For several reasons it may be advantageous to introduce the coordinate q(t) — Q(t) — Q m for 
which the local equation of motion writes 

Mq(t) + iq (t) + Cq(t) = -q cxt (t) (13) 

This transformation is readily performed considering (|7|), with the Qo being replaced by Q m , 
of course. It implies to work with the gg-response function 

(14) 

<Jcxt=0 

instead of the Xcoii(^), with a similar relation leading from the Xosci 1 ^) to the oscillator response 
Xosc(^) for q-motion. The latter is then specified by the transport coefficients 7 = k~ 2 j F , M = 
k~ 2 M F and C = k~ 2 C F . 

Whether or not such an ideal situation is actually met depends largely on the value of 
the coupling constant k. In fair approximation one may say it to be given whenever in (|i~0|) 
the static stiffness C(0) is small compared to the static response x(0)- I 11 such a case the 
strength of all the solutions of the secular equation 1/k + x{ u ) — will concentrate in the 
low frequency domain. Therefore it is important to realize that the static energy may depend 
strongly on temperature. For nuclear physics, for instance, the C(0) drops dramatically above 
kftT = 2 MeV where the influence of shell structure has diminished (see e.g.||25|| ). 

In the more general case the strength distribution of the collective modes, which is deter- 
mined by the imaginary part Xq q {ui) °f the gg-response function, will spread over several (or 
even many) modes. In this case the equation of motion for g(t), which may be obtained from 

(Xqq(u)y 1 q{u) = -g«t(w) (15) 

after a Fourier transform back to time, will contain terms being non-local in time which 
cannot be expressed simply by a sum of derivatives of q(t) up to second order. If one still 
wants to stick to differential form one needs to apply a further approximation. For instance, 
one may simply expand the (x g g(co , )) _1 in ( |15|) to second order in u [^]. Another possibility is 
found in "reducing" the Xgqi 00 ) f° the Xosc(^)- By this we mean to replace the former by the 
latter in a certain range of frequencies. In practice this may be done by fitting the imaginary 
part of Xosc(^) to that peak in x'q q (u}) the corresponding mode of which one wants to treat. 
Notice please that both of these approximations require that only frequencies are relevant 
which lie in some limited regime. In this sense the second one may be considered superior to 
the former. In the past, various computations of transport coefficients have been performed 
which lead to values being in fairly good agreement with experimental experience, see [fZ3]], 



^XcollM = Xqq(u) 



5q(u) 
5{-q ex t{u) 



2.3 Unstable modes 



The oscillator response is analytic in the stiffness C. As we will see below, the latter becomes 
negative whenever the static stiffness C(0) becomes negative, which happens in the region of 
potential barriers where the system is unstable, in the sense that one of the two fundamental 
solutions for average motion begins to grow exponentially. Such a situation can still be 
handled within the locally harmonic approximation as long as the forces do not change too 
rapidly with the coordinate. Recall please that at some given Q the harmonic solution is 
meant to portray the true one only within a limited time interval St anyhow. This growth of 
the time dependent solution reflects itself in the feature that for negative C one of the poles 
of the collective response moves into the upper part of the frequency plane. (Notice please 
that the instability we have in mind still leaves both inertia as well as friction positive). This 
does not cause any problems when dealing with the functional form of the oscillator response 
as given by flT2D, but it may well be so for more general expressions. For instance, one may 
want to write the response function in terms of integral representations, like 

Xixviv) = / — forlmu;>0. (16) 

In this case some caution is required in defining the contour C. For stable systems it is taken 
along the real axes with the poles lying below. This choice is required by causality. If the same 
principle is still to be valid for the case of an instability the contour C must still be chosen to 
be above the one singular pole mentioned before. Actually, besides this "retarded" function it 
may become desirable to have an "advanced" response function as well, for which the contour 
then has to be modified to lie below the poles which correspond to motion backward in time. 
More details of these formal techniques are described in Appendix A. 



2.4 Implications from sum rules 



The integral representation (|T6| ) may serve as a starting point for deriving and exploiting sum 
rules, even in the case of instabilities. Generalizing the common definition to 

s $ = L > ( i7 ) 



the same relations are recovered for the unstable case as those known for the stable one, if 
one only observes the appropriate choice of the contour. For example, for n = — 1 the value 
of the integral in (|T7|) gives the static response, Sj^ 1 ^ = x^v(^ = 0), while for n = 1 one may 
employ the usual arguments to obtain 

S« = I "J*t(n> = Jim - J f = - Jim ^ M . (18) 

One way of getting the equality in the middle is to involve the geometrical series for the 
1/(1 — (fl — ie)/u), which is possible whenever the Sj%) remain finite for n > 0. The latter 
restriction is not really necessary. One may simply replace the Q in the first integral by 
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Q = uu(l + VL/uj) = [1 — (Q/uj) — (Q/uj) 2 / (1 + Q/uj)]^ 1 and observe that for uj — > oo the last 
term in the brackets gets less and less important. Then one only needs to have the Sffj vanish, 
which in our case is no problem as later on we are going to apply ( |T8| ) only for the diagonal 
case. Please observe that for the oscillator response ([12]) the validity of the expression on the 
very right follows directly from inspection. 

The sum rules for n = =)=1 are intimately related to two of the three transport coefficients. 
This can be inferred from the values these two sums take on for the oscillator response, namely 

SL' ] = ^ S oll = ^ with 5ffl = 0. (19) 

For later purpose we may note here that these results remain unchanged when the oscillator 
is modified by introducing a frequency dependent friction 7 = 7(0*). 

As the oscillator response is derived from the original, microscopic expression given by 
(H) together with (|14D, the sum rules may be used to derive inequalities for the transport 
coefficients. How these inequalities come about may best be seen writing the full response 
function in terms of a pole expansion like 

V 2M a £ a \uj-uj+ uj-uj a J 

with = ±£ a — iT a /2. Then one may deduce results like those given in ( |T9| ) for each term 
of the sum over a. It is then evident that the value of the S^ 1 ^ for each term cannot be larger 
than the corresponding value of the total sum. 

Let us look at the inertia first. The form (^) implies lim^oo w 2 (x co \\{^) — x(^)) = 0. 
Therefore we get 

k- 2 S%= ^yUm= ^x"mn = U[F,F]) (21) 

J-00 7T J-00 TX II 

and consequently 

^ < S« = ^<[>^]> = - (22) 
M n mo 

which means that the inertia for the one mode we have chosen to pick cannot be smaller than 

the mo given by the microscopic value of the energy weighted sum. In nuclear cases this mo 

can be seen to correspond to the inertia one would get in the liquid drop model for irrotational 

flow (see e.g. 0). 

Next let us look at the stiffness. Here the situation is slightly different as the mode we 
like to look at explicitly may become unstable. However, assuming that all the other ones 
remain stable one may simply repeat the arguments from before such that for S^ 1 ^ = Xgg(0) 
one gets 

1 < y (Q) - k X{0) - 1 X{0) (23) 

These expressions for Xgg(0) are obtained after observing (|]), (|14D and (pTOf ) . The last equation 
shows the Xggify to have the same sign as the static response C(0), from which observation 
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it follows that C will become negative as soon as (7(0) is negative. In these arguments we 
have assumed the static intrinsic response x(0) to be positive, as there is no reason to expect 
the intrinsic system to be unstable. As mentioned before, we want to stick to the case where 
this value is larger than the absolute value of the static stiffness. The inequality allows one 
to prove that effective stiffness cannot be smaller than the static one: 

C > (7(0) (24) 

For stable modes this is immediately clear, for unstable ones it is useful to make a detour and 
introduce the absolute values (7 = — | (7 | and (7(0) = — | (7(0) |, in which way one gets from 



3D 

111 1 

\c] ~ YcW\ i- I c(o) | / x (o) - |cfi ' (25) 

The inequality ([24]) expresses in clear fashion that in general the effective stiffness (7 will be 
different from the static one, (7(0). The reason for this difference is found in dynamical effects 
contributing to the effective potential. It is interesting to note that at an instability these 
effects actually lower the absolute value of the stiffness, opposite to the case in stable modes. 
In this sense the dynamics weakens the instability by de-creasing the conservative force. 

2.5 Inferences for the high temperature limit 

Two body Hamiltonians of the type (^]) are frequently used in nuclear physics. Commonly, 
however, the notion of the coupling constant k is just taken in literal sense, namely to be 
a given parameter which does neither change with deformation Q nor with temperature T 
(see e.g. the review [p9| ). This is an essential difference to the situation described above 



for which k does depend both on Q and T, as given by fllPl). To see the implications of the 
variation with Q would require more extensive studies, which are to be deferred to the future. 
In particular, it would be necessary to extend the present treatment to the multi-dimensional 
case, as given for genuine quadrupole vibrations, for instance. Such a procedure is possible, 
in principle (see e.g. WW), but it involves efforts much beyond the present study. For the 



T dependence, on the other hand it is quite easy to see important consequences, even at a 
qualitative level by simply looking at the limit one would get for the collective mode when 
increasing the nuclear excitation. 

In the last decade or so, rather fancy procedures have been adapted from field theoreti- 
cal methods to study collective dynamics at finite T. Some of them exploit the concept of 
small amplitude vibrations about some given thermal state of the system, like RPA at finite 
temperature (see e.g. [[29]]) or the variants of the "Static Path Approximation" such as the 
P(erturbed)SPA of |30| or the C(orrelated)SPA of ||31|| . Also there separable two body inter- 
actions similar to the one given in (0) are exploited (with the k being some given constant, it is 
understood). One of the striking results is that with increasing T the vibrational modes turn 
into the unperturbed excitations of the nucleonic degrees of freedom, even for those which 
at T = would show the typical collective behavior. In our opinion, this feature contradicts 
basic concepts of Niels Bohr's model of a compound nucleus. The latter is dominated by 
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the effects of strong (incoherent) two body correlations, or, to use a word which draws more 
on analogies with transport theory, by the effects of two body collisions. There can be little 
doubt that the latter will become the more important the larger the thermal excitation (c.f. 
the Ansatz for the nucleon's self-energy suggested by (|TT|)). One may then legimately ask, 
why it should be possible for the independent particle motion configurations to dominate the 
dynamics at large T? 

At this stage it may be very worth while to recall the essence of the Strutinsky renormal- 
ization of the total static energy For small excitations the variation of the latter is completely 
determined by shell effects. However, what remains when they are gone is the liquid drop 
energy. Consequently, the vibrations of the nucleus should behave more like those of a " drop" 
of nuclear matter showing largely "macroscopic" behavior, rather than that of independent 
particle motion. This is indeed what one gets from the theory described in the present paper 
(see e.g.[[2lj or [12|], with further references given therein). As mentioned previously, with 
decreasing shell effects the coupling constant drops to its macroscopic limit. Consequently, 
the modes get very soft. Even more important is a strong concentration of strength in just 
one broad collective mode. The large width associated with this has to be understood as 
clear signature of strong damping, the latter feature being inherently related to the strong 
correlations mentioned previously. Unfortunately, the strength distribution of isoscalar modes 
has not yet been measured for finite temperature. Such a measurement would be an ideal 



test of the transition to the macroscopic limit [25], |12|. However, this transition reflects 



itself in transport coefficients for slow large scale collective motion like fission, in particular 
in their temperature dependence. Within the theory described above one finds that damping 
increases with T, which is in at least qualitative agreement with experimental findings, see 
[ f2"5|] and \T1\. It should be noted that the relevant quantity is not the friction coefficient 



itself, but combinations with those of stiffness and inertia, depending on which quantities one 
looks at. Whereas the T-dependence of the stiffness is well known and can be deduced from 
calculations of the static energy, this is less clear for the inertia. It has been shown first in ]32 



that one of the consequences of the T-dependent coupling constant k of ([TOj) is a behavior of 
M(T) as one would expect for a transition to the macroscopic limit: At small excitations the 
M shows the typical features of the cranking model, but above the point where shell effects 
have disappeared M approaches values similar to those of irrotational flow. 



3 Extension to dynamical fluctuations 

In the more general sense our treatment so far has not gone beyond a mean field approxi- 
mation. It is true, however, that the dynamics in the basic quantities (F) t or q(t) may be 
damped. Microscopically this damping implies that we need at least to consider effects which 
are not contained in a pure Hartree (or Hartree-Fock) approximation. However, it still suffices 



to stick to a description in terms of self-energies in the single particle degrees of freedom |12 
and the q(t) still simply is a c-number. Evidently the latter restriction has to be given up 
whenever one is interested to describe processes where fluctuations in the collective degrees 
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of freedom become relevant. Take the example mentioned before, the decay out of a potential 
minimum. One knows from the pioneering work of Kramers |J how in the realm of classical 
physics the dynamics of fluctuations across the barrier can be treated and how they influ- 
ence the decay rate. In |33|] it was shown how the quantal extension of this decay rate can 
be obtained within a description based on transport equations. Some of the ideas behind 



this derivation were borrowed from the method developed in [18]. Unfortunately, also this 
method was again based on functional integrals for which reason it can be applied only to the 
schematic cases mentioned in the introduction. In this section we thus would like to proceed 
describing the basic steps needed for a formulation in terms of transport equations, with an 
elaborate discussion to come in the last two sections of the paper about the essential ingre- 
diences of this transport equation. By those we mean these quantities which are intimately 
associated to fluctuations, namely the coefficients which parameterize the diffusion process in 
collective phase space. 

3.1 Quantization of the collective variables 

As the first step we need to render the collective degree of freedom to become a fluctuating 
variable. Often people are tempted to do this simply by adding to the Hamiltonian (^J) a 
kinetic energy term like IT 2 /(2m), with some unknown mass parameter m to represent the 
collective inertia and II being the momentum conjugate to the coordinate, which then is to 
be considered an operator Q. It is not difficult to convince oneself that in this way self- 
consistency would badly be lost; for a detailed discussion see fl2]] . This feature cannot be 
cured by playing with the potential part, which means to say by adding some function of 
Q — Q Q . Rather one needs to introduce a momentum dependent coupling. 

Starting from the effective Hamiltonian (|2|) in the original particle space, this method 
allows one to deduce the following Hamiltonian for the total system (see ||11||): 



H = H(xi,pi, Qo) +kUF-^(Q- Qo)F + + ^-^-(Q - Qo) 2 (26) 

Comparing with ([]]) one easily identifies the parts for bare intrinsic and collective motion, as 
well as the coupling between both, which here consists of the two terms being linear in n 
and Q — Qq. In (|26"[) the F = i [H(xi,pi, Qo), F], the n stands for the canonical momentum 
and mo is the inertia which by (p2|) was introduced to represent the energy weighted sum, 
and which serves here to define the unperturbed collective kinetic energy. The parameter (3 
has to be specified only in connection to the collective momentum, a point to which we will 
come to below. This f3 drops out of any equation of motion which only involves the collective 
coordinate Q. This feature is very similar to the properties we found before for average 
dynamics (as expressed in (F) or in q) which could be traced back to the response function 
([3]). It should be no surprise that the same function appears again when the dynamics for 
average motion in q(t) is re-derived from Ehrenfest's equations to (|2~6D . We may note in 
passing that in the original version of the Bohm-Pines procedure, as it had been developed 
for the electron gas in 1(J , no such additional parameter /3 was needed. The difference to our 
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case lies in the fact that we want to apply the method to an attractive two body interaction 



in contrast to the repulsive one of the electron gas (c.f. [11 1). 

The description in terms of the Hamiltonian ( p6|) is not yet complete. As it stands the 
H would have too many degrees of freedom and in this sense could not be equivalent to the 
original two body Hamiltonian of (0). Indeed, the Bohm-Pines procedure necessarily 
involves a subsidiary condition which for the present case reads 

kF = Q-Q (27) 

It is nothing else but the "operator version" of the self-consistency condition (]?]). This property 
and the very fact that it is inherently incorporated into the form (|26|) of the Hamiltonian H 
makes it very plausible that the description of collective motion on average is the same as 
before. The new feature is seen in the fact that now we are in a position to treat fluctuations 
in the q in the way similar to what the Hamiltonian of (|2[) allows one to do for the 
microscopic quantity F . This feature goes along with the fact that in fl2"E| ) no two body 
interaction appears anymore. Rather, its effects are hidden in the terms involving collective 
coordinate and momentum. 

With the Hamiltonian (|26|) one may proceed to derive effective equations of motion for 
the propagation in collective phase space. In this spirit the Nakajima-Zwanzig projection 
technique has been applied in |7j to get the transport equation for the Wigner function of the 
collective density. However, to account for self-consistency not only on the level of average 
dynamics, or on that of the mean field, a non-perturbative method had to be developed for 
treating the integral kernel. Again this was possible within the locally harmonic approxi- 
mation for which it suffices to know the propagators for the consecutive time steps St with 
St <C r co ii. As for St — > this propagator starts at a definite point in phase space, it was 
argued that a Gaussian approximation would do. Thus one only needs to specify the time 
evolution of the first and second moments. For this non-perturbative Nakajima-Zwanzig equa- 
tion it was possible to prove that it is in accord with the (quantal) fluctuation dissipation 
theorem (FDT) 0], which in full glory requires to have non-Markovian equations of motion. 
This feature makes it possible to take advantage of the FDT for establishing the approximate 
equations of motion, which for the first and second moments may be taken to be of differential 
form. 



3.2 Harmonic approximation for first and second moments 

For average dynamics we have seen the differential equation of motion to come up through the 
reduction of the "full" response function (0) (or the one of (0)) to that of the oscillator given 
in (fT2|). Evidently, one would like to have the equations of motion for the second moments 
to be of similar type. Clearly, in addition to the coordinate they will involve the collective 
velocity or the corresponding momentum as well. Let us first introduce the latter on the level 
of first moments, i.e. for q c (t) = (q) t and p c (t) = (p) t - It is clear that for the harmonic case 
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we are looking at, one expects the following set: 

M J t ^(t) = Pc(t) j t p c (t) = ~lj t <lc{t) - Cq c (t) , (28) 

As the first equation tells one, the momentum appearing here is the kinetic one, called p 
henceforth. It is here where one may benefit from having the parameter {3 appearing in (|26|) . 
As shown in [|TTJ and JFj, the choice f3+k = k 2 C allows one to replace the canonical momentum 



IT by the kinetic one not only for average motion but also for the fluctuations around it. As 
demonstrated in J7] these fluctuations satisfy the following linear set of equations 



-S (t)-* 
dt qq{t) M 



E qg (t) - ^ qp {t) = (29) 



I 1 

-X qp (t) - jjX pp (t) + CT, qq {t) + jjZ qp (t) = D qp (30) 



-E pp (t) + 2CY, qp {t) + ^E pp (t) = 2D PP (31) 



— Spp(t) + 2CL qp {t) + 

The second moments are defined as Tl qq (t) = {q 2 ) t — Q 2 (t), and similarly for the combinations 
{q,p} and {p, p}. It is seen that the structure of the homogeneous parts follow from that 
of for which reason the q 2 {t) etc. satisfy the set ( P9"Q to ( |3T| ) with the inhomogeneities 
D qp , D pp put equal to zero. For any finite value of these diffusion coefficients D qp and D pp the 
fluctuations will become finite in the course of time. This will be so even for zero initial values, 
as it is the case when the solutions are to represent the second moments of the propagators 
(see above and below). So far the diffusion coefficients have not been specified. It is here 
where we may benefit from the FDT. As is easily verified, these diffusion coefficients can be 
expressed by the stationary solutions of (p9|) to ([31]). Let us first take the conventional case 
of a positive C for which these stationary solutions are those of equilibrium such that we may 
write: 

p ~ M pp ~ OZj ot M p p ( ' 

In the next chapter we will exploit the (quantal) FDT to evaluate these expressions. After- 
wards we will turn to the case of an instability with C < which will be handled by analytic 
continuation. 

Often one encounters cases of strongly over-damped motion. It is not difficult to convince 
one-selves (see ]12| and c.f. 0) that in this case the set of equations reduces to (|29| ) to (|3T| ) 



jfM + ^ qq (t) = 2DZ d (33) 

with the diffusion coefficient D° v q d = 

It may be helpful to the reader to add some further comments on the seemingly Markovian 
nature of our basic equations of motion. There can be no doubt that approximations of this 
kind are necessary in one way or other. However, for average dynamics this step to differential 
form does involve less stringent conditions than one would expect on usual grounds. This 
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feature is intimately related to the construction of the oscillator response or the definition 
of the transport coefficients M, 7 and C. Look once more at (0) and suppose that the 
microscopic response x(°°) can be represented by a Lorentzian, or its dissipative part, rather, 
which often may portray a realistic situation quite well. For such a case the form (O) follows 
without any additional approximation. On the other hand, it is almost evident that the set 
of equations ( P§| ) to ([H]) for the second moments can hardly be correct for all times t and 
in particular not for very small t. This set supposes that we use it for large times, a feature 
already suggested by our way of calculating the diffusion coefficients: They are chosen in such 
a way as to warrant that the set ( f29|) to ( j31~l) describes correctly the relaxation to the proper 
equilibrium. 

3.3 Treatment of large scale motion 

The information contained in fl28|) to ([H]) may now be used to construct the time evolution 
in global sense. Two possibilities come to one's mind. First of all and as indicated before, one 
may take the first and second moments to define a Gaussian. The latter may be constructed 
such as to represent motion within the time lapse 5t = £2 — £1 and which at t = t\ starts from 
some point Qi, Pi in phase space. For a sufficiently short interval St the widths T.^t) may 
be assumed to remain small such that this motion remains restricted to a narrow region in 
phase space, thus justifying the linearization procedure we started off. As this propagation 
can be performed for any point Qi, P\, one may in the end sum up all possibilities needed 
to construct the time evolution of any initial distribution and even over time intervals of 
macroscopically large size. That this "propagator method" may be used even to describe 
such complex problems as the decay out of a potential minimum has first been demonstrated 
in @. 

Numerically this method turns out to be time consuming as the "summation over all 
possibilities" requires integration over the phase space not only for the first step but for each 
following intermediate step as well. It so turns out that it is simpler to translate it to the 
picture of a Langevin force. There one follows individual trajectories and constructs the 
final distribution by repeating this latter procedure for a suitably large manifold of points 
representing the initial distribution (in global sense). These trajectories are constructed in 
such a way that for the elementary time step 5t = t2 — t\ one first moves along the average 
trajectory, say by solving the set (p8|), to find the point representing the actual position at 
time £2 by Monte Carlo techniques. This sampling is to be performed by taking into account 
the size the fluctuation attains at t = £2 according to the fluctuating force, which in turn may 
be determined from the information contained in the diffusion coefficient. 

4 Diffusion coefficients for stable modes 

To evaluate the diffusion coefficients from we need to calculate the values the fluctuations 
would take on if the (collective) system were in equilibrium locally. The fluctuation dissipation 
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theorem implies 

S-=^/^coth^^». (34) 

Here the indices u, \i are supposed to represent q and p. Exploiting the relations Xqpi^) = 
iMu>x'qq(u) and Xp P (uj) = (Mu) 2 x'q q (uj), all fluctuations can be expressed by the one response 
function Xqq- The different behavior of q and p under time reversal causes the cross term to 
vanish, i.e. E^ = 0. Formally this is reflected by x'q p being even in uj. It is for this reason, 



that no diffusion coefficient appears in ( 29|) 



Actually, the values given by (|34|) are not yet suitable for our purpose. It was mentioned 
before that the equations ( p9|) to ( |3~1~D should be consistent with those for average motion 
given in (p8|). As the latter are derived from the oscillator response we are forced to take the 
latter also when evaluating the fluctuations from an expression like (|34]). In other words, the 
diffusion coefficients ought to be calculated from relations like fl32|) but where the E^, are 
replaced by those for the oscillator E°^ c , 

7 osc osc 1 

D PP = jTf^PP D QP = C^° qq C - J^^PP ' ( 35 ) 

and the equilibrium fluctuations are finally to be calculated from 

S - = h J d ^ C ° th {^) X ' L{uj) (36) 

du I hiv\ 
— coth — 
2tt \2T J 

In doing so we keep the E^ = untouched, as it results from a fundamental symmetry. 



and 

-iosc _ )\#2* / ""^ / ' tu/ \ , .2 ,// 



Kp C = M h l ^ coth I 2T J u ^osc(^) • (37) 



4.1 Limiting cases 

Before turning to the general case let us look at the limits of high temperature and of small 
damping. The former is established if the argument of the hyperbolic cotangent is small, 
i.e. hu/(2T) « 1, for all frequencies which contribute to the integrals. Quite generally, 
one finds: E^ = Tx vt i{w = 0) which for the oscillator leads to the classical equipartition 
theorem: 

( P 2 )eq _ T _ C( Cf > T (3g) 



2M 2 2 

As a consequence, for the diffusion coefficient one obtains the classic Einstein relation 

D pp = 7 T with D qp = (39) 

Notice that for this result no assumption was necessary about the strength of the friction 
coefficient, besides the fact that 7 will influence the range of frequencies which must be 
considered small compared to T. 
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In the limit of zero damping one finds a similar result, only that the temperature T has 
to be replaced by an effective one, namely 

T*(^) = — coth(^— J (40) 

with to 2 = \C\/M. This is easily verified letting the friction 7 in the oscillator response (O) 
go to zero. Then the dissipative part of the response function can be written as 

hm = (% - w) - % + t*)) (41) 

The resulting version of the generalized Einstein relation, i.e. D pp = 7T*, has been used 



pp 

in 



26| . It allows for a nice and simple evaluation of the correction terms to the classical 
relation, namely D pp ph jT(l + (l/3)(ctj/2T) 2 . However, as we will demonstrate below, the 
form D pp = 7T* approximates the correct formula only for extremely small friction. 

4.2 General case 

To evaluate the integral in ([34]) one commonly expands the hyperbolic cotangent into the 
uniformly convergent series 

, (Tiu)\ 2T 1 ^ 2ttT , , 

coth — = — > = — — 42 

\2T) h ^w-ine h V } 

and applies a general relation (see (|A.5|) ) which allows one to evaluate the response function 
through its spectral density, the dissipative part; such expressions are briefly discussed in 
Appendix A. For the terms with positive n, the integral leads to the retarded response function 
calculated at the corresponding pole of the coth, namely 



du = R 

c 71 u — zlnlG ' '" 



X5,(i|n|e) (43) 



For negative n, one may use a relation between retarded and advanced response to find 



c 7r uj + i\n\Q 



(xi(i|n|6)) (44) 



Notice please that in the classical limit it is only the term with n = which survives and 
which leads to Tx vii {uj = 0). The quantum corrections are caused by those contributions to 
( f42| ) having a finite n. 

A straightforward calculation then shows the diagonal case elements to be given by 

m = T ( Xlllt («> = 0) + 2 £ Xw(ine)\ (45) 

In this expression we have used the fact that the response function \nn is re& l along the 
imaginary axis. 
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To evaluate the diffusion coefficients we need to apply ( fl5|) to the oscillator response. 
Unfortunately, the resulting series diverges logarithmically for the momentum fluctuations. 
Actually, such a behavior can better be read off from the integral fl57|) after recognizing that 
the dissipative part of (^) decreases with u> only like u~ 3 . To cure this problem one needs to 
regularize the integral. In J7| the latter was simply cut off at a certain upper limit. Here we 
would like to apply the somewhat more elegant method of the so called Drude regularization 
(see e.g. ||). It implies replacing in (|l|) the friction coefficient by a function of the form 

^ = Y^k^ ■ (46) 

Properties of the resulting modified oscillator response function are discussed in Appendix 
|B[ In ( f46|) the additional parameter wd appears, the appropriate choice of which depends on 



general properties of the system. In the present work we shall assume the time scale 1/vjd 
to be well separated from that of collective motion and will examine how in a physically 
reasonable regime this parameter influences the diffusion coefficients. 
For the £°p C one finds the following expression: 

S^ c (w D ) = TM + 2 fj [C + nQj(inQ)] X o SC (inQ)j , (47) 

A common procedure to evaluate further both (|45|) and (|47|) is to make use of the digamma 
function ^(z) = ^lnT(z). It allows one to sum up series over products of pole terms, for 



instance in the form (see e.g. |27, 2£] ) 



y 1 = tt(l+y)-*(l + z) (45 

il (n + y){n + z) y-z 



The oscillator response with the constant 7 replaced by the one of ( [46]) has three poles 
(see Appendix [B]). Thus one may rewrite it as a sum of three terms where each one factorizes 
like the left hand side of fl48|). After some manipulations one then ends up with 



c K ^ D) = r (i + g * ( 1 + ^ I ) m 




and 



M pp {Wd) = °^ {Wd) + ~M7 g M,)-^))(^)_^ +2 )) *{ 1 + l 2^) (50) 



where the convention 

^(i+3) = is 

used. In ( ^9|) also the gg-fluctuation has been evaluated 
with the Drude regularization, mainly to allow for a clear comparison to the other case where 
this regularization is really necessary. It is clear from (HB) that for large wn the result (PE3) 



will turn into the one calculated for constant friction. How rapidly this happens depends 
on the values of the parameters involved. This is demonstrated in Fig.l a by numerical 
calculation. It is seen that large values of vje> are needed for small values of temperature and 
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Figure 1: Equilibrium fluctuations in q (upper panel) and diagonal diffusion coefficient 
(lower panel) as functions of the cut off frequency for various values of damping rate r\ and 
temperature T. (The calculations are performed in units with k B = 1 and Ti = 1). 

large damping before the asymptotic limit is really reached. It may be seen as well that even 
for comparatively small zud the difference to this asymptotic limit is quite small. (It should 
be clear from above that wb must not be chosen smaller than the typical frequencies of the 
original oscillator.) Here and in the following it turns out convenient to measure frequency in 
units of w = \J\C\/M and the degree of damping by the dimension- less ratio rj — 7/ (2Mzu). 

Different to the gg-fluctuations, those of the momentum would diverge for large zud- 
Therefore, one must expect them to exhibit a stronger dependence on the Drude frequency 
zu D than seen in the qq— fluctuations. This is demonstrated in Fig.l.b. Indeed, for large 
rj this dependence is quite sensible. We may notice that the diffusion coefficient increases 
with zud, although its classical limit is independent of zud- In this sense it can be said the 
quantum corrections to increase with zud, too. The standard values chosen for the figures 
presented here, namely vdd 10 u, are probably very reasonable. Fortunately, in this region 
the dependence of the diffusion coefficient on wu is not very strong, the value of D pp changing 
by 30% if wd is increased by a factor 2. 
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Let us turn now to the dependence of the diffusion coefficients on temperature T and on 
the effective damping rate rj. Although the transport coefficients for average motion, M, 7 
and C, must be expected to vary with these parameters as well, this feature will be neglected 
here. All computations will be done in Drude regularization. If not stated otherwise, the 
value of zud will be chosen such that the ratio ( = vj/wd equals 0.1. To simplify notation 
we want to use the same symbol as if we would use the mere oscillator, which is to say that 
in expressions like the ones on the left hand side of (f49[) and ( |50|) the argument vj d will be 
left out. We shall also omit the index osc at times. 

Let us begin looking at the equilibrium fluctuations (see also f|), first in the coordinate 
q. Rather than showing the T,° q c itself, we plot CH° s q c , normalized to its value as given by the 
quantal ground state (Fig.2.a). Fig. 2. a shows the temperature dependence for different values 
of rj. One finds the quantum effects to disappear both with increasing T and rj. The first effect 
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Figure 2: Equilibrium fluctuations in the coordinate (upper panel) and in the momentum 
(lower panel) as functions of temperature normalized to the frequency w of the oscillator, for 
different values of the damping rate 77. 

is expected, although it may be somewhat of a surprise how quickly the high temperature 
limit is reached. The second effect causes the to reach this limit the faster the larger the 
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damping. While behaving similarly with increasing T, the fluctuations in momentum differ 
in their dependence on rj: They increase with increasing damping. This is demonstrated 
in Fig.2.b. Here £°p C /0.5Mm (or equivalently D pp /0.5-ym) is plotted versus temperature for 
different rj. This figure clearly shows how much the quantum effects are enhanced by friction. 
Take the regime of T < m where the change with T is quite small. For large r\ the fluctuations 
exceed several times the values they have for undamped motion. 

Next we study the diffusion coefficients. They are plotted in Fig. 3 by the dashed lines 
as function of r\ for different temperatures, the latter being normalized to m. In Fig. 3. a we 
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Figure 3: Diagonal (upper panel) and off diagonal (lower panel) diffusion coefficient for 
stable (dashed) and unstable (fully drawn) modes, as functions of the damping rate rj, for 
three different temperatures. In the lower panel the curve for T = 0.1m shows 0.1D qp /T, as 
indicated by the symbol xO.l. 

show the diagonal one, divided by its high temperature limit, D pp /jT; Fig.3.b represents the 
non-diagonal one, divided by temperature. Both figures demonstrate that for the diffusion 
coefficients quantum effects increase with damping, and that they may become quite large. 
Take the non-diagonal element plotted in Fig.3.b and recall its definition according to (|35|). 
It measures the difference between CE° s q c and jjT,° s p , which vanishes in the high temperature 
limit. But Fig.3.b demonstrates that for small T this difference, and thus D qp , reaches values 
which are several times this unit. Moreover, these figures show that the "limit of zero friction" 
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discussed previously must be taken literally: It is obtained for such small values of friction 
only that it is practically never realized in nuclear physics. Look for instance at Fig.3.b: 
Already for rj > 0.5 the D qp is seen to be at least of order T even for T « w, whereas in the 
strict limit of zero friction D qp vanishes identically. 

It has been mentioned that in the computations of ]7J regularization of the appropriate 
integral had been achieved by introducing an upper limit in the integral itself, rather than 
to change the integrand applying the Drude regularization. The two methods are compared 
with each other in Fig.4. They are seen to lead to the same dependence on 77 but may require 
an optimization of the cut-off parameters, if the comparison is to become more quantitative. 
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Figure 4: Comparison of two regularization schemes: Sharp cut off (thin line) and smooth 
cut off through Drude regularization. The values of the cut off frequencies, w max and u D 
respectively, were choosen in such a way that in both cases the diffusion coefficient attains 
about the same value. 



5 Diffusion coefficients for unstable modes 

The principle which allows one to deduce diffusion coefficients from the fluctuation dissipation 



theorem even for unstable modes has been suggested in ||, see also fl55| . In the following we 
like to explore this possibility in more detail with the help of formulas obtained by adapting 
linear response theory to instabilities. All we need to do is to generalize expressions such as 
( [45D and fl5H| ) to the case of negative stiffness. 

The basic idea is found in realizing that all the expressions derived and discussed in the 
previous sections are analytic expressions in the transport coefficients, and thus also in the 
stiffness coefficient C. This is true in particular for the general expressions for the diffusion 
coefficients as given by (f49"D and (|50|) , together with fl35"l) . In Fig. 5 we present a computation 
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of three typical quantities, namely CE^, D pp and D qp , all three normalized to temperature, as 
function of the stiffness parameter C (divided by the inertia M and relative to (7/M) 2 ). The 
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Figure 5: Diffusion coefficients (diagonal: thick fully drawn line, off diagonal: dashed line) 
and analytic continuation of the equilibrium fluctuations in q (thin full line) as functions of 
the stiffness C, for fixed values of M and 7. 



result demonstrates the smoothness, or analyticity, across the point C = 0, which actually 
corresponds to the inflection point C(0) = of the static energy. This feature is quite 
interesting, for several reasons. First of all one realizes that the fluctuation in the coordinate 
diverges when C approaches zero from above: X^j — > +00 for C — > + . Secondly, across this 
point this "fluctuation" gets even negative. Fortunately, for the equations of motion it is not 
this " hypothetical" static fluctuation which matters but the diffusion coefficients. The figure 
nicely demonstrates that all of them remain well behaved at this point implying that the same 
will be true for the dynamical fluctuations! 

On the other hand, Fig. 5 also demonstrates that there are limitations in applying the 
analytical continuation to get meaningful diffusion coefficients. It should be clear that the 
diagonal one ought to be positive. However, depending on the parameters involved, like 
the transport coefficients for average motion or temperature, the D pp may become zero or 
even negative. This feature can be seen more clearly in Fig.6 where we plot the analytical 
continuations for and as function of temperature divided by the bare frequency, i.e. 
of T/w. The two quantities are normalized such that in the high temperature limit they 
approach the same value, namely 2T/zu. Actually, the factor two here is a relict of the 
discussion of the stable case presented in Fig.2, where the same quantities were plotted, and 
where this normalization was convenient at T = 0. Since the curves of Fig.6 just represent the 
analytical continuation of those in Fig.2, it is not surprising that the quantum effects behave 
in the same way as for stable modes when studied as functions of temperature and friction. 
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T/G3 

Figure 6: Analytic continuations of the (diagonal) equilibrium fluctuations to the case of 
negative stiffness (dashed line for the coordinate, fully drawn lines for the momentum), as 
functions of temperature for different values of friction. The insert shows the momentum 
fluctuations for a larger range of the temperature and the same values of r\. 

At this stage it may be worth while to refer the reader to the two easy examples discussed 
in the previous chapter, namely the limit of high temperature on one hand and the case of 
zero damping on the other one. Indeed, eq.(|39|) tells one that the relevant diffusion coefficient 
D pp remains positive for all T, but from eq.(p8|) on sees the to change sign at the point of 
inflection. On the other hand, the T* gets negative below a critical temperature of T c = hw/it. 
The figure Fig. 6 shows such features to hold true in the general case with the following two 
important modifications: a) The T c is shifted to lower values when rj is increased; b) this shift 
is bigger for the momentum fluctuations. At larger values of 77 the latter even show a different 
behavior, in the sense that they become larger than the value given by the high temperature 
limit, which for larger T they approach from above. 

The non-diagonal diffusion coefficient D qp is plotted in Fig. 7 as function of temperature, 
both for positive and negative stiffnesses, and in Fig.3.b by the fully drawn lines as function 
of friction. The values of D qp for positive and negative stiffness differ from each other only 
for small temperatures and small friction, and practically agree with each other already for 
temperatures larger than about 0.5hw and friction 77 > 2. A similar behavior is true for the 
diagonal element, as may be inferred comparing fig. 3. a and Fig. 6. At very small T jw the D qp 
shows an anomalous decrease with Tjw (see Fig.7 ) or with 77 (see Fig.3.b ). This actually 
happens in a regime below T c , the reason of which will be discussed on below. 
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Figure 7: Off diagonal diffusion coefficient as function of temperature for three values of 
friction, for a stable (dashed lines) and an unstable (fully drawn lines) system. The curves 
for r] = 1, 10 show D qp multiplied by 0.2 and 0.1, respectively. 



Finally, we would like to come back once more to the fluctuation dissipation theorem 
in its general form. We recall that it was this theorem which allowed us to determine the 
diffusion coefficients in terms of the response function of average motion. Furthermore, these 
are exactly the expressions, taken as analytical functions of the stiffness C, which we were 
able to generalize to unstable modes by analytical continuation from a positive stiffness to a 
negative one. The construction of the special contour C now allows us to perform these steps 
in a more direct fashion, provided one more constraint on the contour is considered. That 
is to say that at instabilities we may express the "equilibrium fluctuations" by the following 
integral: 

X'M (51) 



yeq 





f CO \ 


[ — coth 




lc 2vr 


v2T/ 



In this integrand there appears another factor having poles in the complex plane, the coth 
which diverges at the Matsubara frequencies u;^ = ±n2irT/h. As in the stable case, the 
contour C has to cross the imaginary axis in between those which lie closest to the real axis 
[ ftB| , namely ±2ttT '/%. This condition actually puts a lower limit on the range of temperatures 
for which such a construction may work, namely T > Tq = K\u+\/(2tt), with u + being the 
frequency of the unstable mode (see Appendix B, in particular eqs.( |B.9| ) and (|B.11|) and 
the text below them). This follows from the constraint discussed in the context of eq. (TO) 
according to which the contour C has to lie above the unstable pole. Actually, this T coincides 
with the so called "cross over temperature" encountered in dissipative tunneling (see eq.(3.7) 
of @), which will be commented on in the next chapter. For illustrative purpose we may 
look once more at the limiting case of zero friction for which one gets T = Tiw /(2tt) and thus 
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2 To = T c . The feature of the critical temperature T c being larger than the To is generally 
correct. 

It is interesting to note that such a minimal temperature below which motion can no longer 
be treated in a harmonic approximation is found also in two extensions of the SPA, called 
P(erturbed)SPA in [3C] and C(orrelated)SPA in PT . Both approaches are very similar to one 
another, in the sense that small amplitude vibrations are considered on top of or around the 
"static path" approximation to the functional integral — and are thus analogous in spirit to 
our locally harmonic approximation. Being based on a path integral formulation, it may be 
not very surprising that in this method the minimal temperature turns out to be identical to 
T . It must be said, however, that both in PSPA as well as in CSPA no damping of the modes 
is considered. Within our time dependent formulation the latter feature is made possible 
through considerations of coupling of lp-lh-correlations to more complicated ones. In this 
sense the two body interaction effectively goes beyond the separable form given in (|2|), which 
just refers to the one generating collective motion. Certainly, more work will be needed to 
make the correspondence of these different procedures more apparent. 



6 Summary and Conclusions 

In this paper we have reviewed a method which allows one to describe the dynamics of 
damped motion in collective phase space. The main emphasize was put on the evaluation of 
diffusion coefficients which govern motion of the second moments both for stable as well 
as for unstable modes. The latter case was handled by a suitable analytic continuation of the 
(quantal) fluctuation dissipation theorem. To this end some elementary properties of linear 
response functions had to be generalized, too. 

From the numerical results it became evident that, when applied to the diffusion coeffi- 
cients, this continuation looses its physical meaning below a critical temperature T c below 
which even the diagonal element D pp would become negative. On the other hand, for the 
fluctuation dissipation theorem, as well as for some basic relations of response theory, such 
a continuation is possible even for temperatures down to a T < T c . This T was seen to 
be identical to T = h\u + \/ (2tt) where u + = — (oo + )* measures the position of the pole along 
the positive imaginary axis which belongs to the unstable mode. This To is thus identical to 
the so called " cross over temperature" Q one encounters in the phenomenon of " dissipative 
tunneling" and for undamped motion in the PSPA and the CSPA of |30| and | 3l| , respectively. 



This association to tunneling is by far not accidental. As mentioned above, the equations 
of motion associated to the harmonic case can be exploited for describing global motion within 
a locally harmonic approximation. In the end this means to describe collective dynamics by 
way of a transport equations like the one of Kramers, with the only difference of generalizing 
the classical Einstein relation to the proper quantal form. This implies to have two diffusion 
coefficients D qp and D pp instead of the single one D pp — jT. Applying this scheme to 
evaluate the rate of decay out of a potential minimum one finds |33| the formula R = /q-Rk 
known from such model cases for which an application of functional integrals becomes feasible. 
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Compared to this latter situation, the derivation sketched here has both an advantage as 
well as a disadvantage. The advantage can be found in the fact that this locally harmonic 
approximation can be applied to cases for which the coupling between the collective degree 
of freedom and the rest is too big to allow for the simplification of linearizing in the "bath 
variables", with the latter simply being represented by oscillators with fixed inertias and 
stiffnesses. On the other hand the description in terms of transport equations of the type 
mentioned is valid only for T > T c > T . In this sense this critical temperature T c indicates 
a breakdown of the locally harmonic approximation in general. 

In this paper it has been emphasized that the construction of the basic equations of 
motion is possible in self-consistent manor. By this notion we mean that the collective degree 
of freedom is introduced to represent the dynamical microscopic quantity one wants to study. 
This is necessary for all systems which per se exist of identical particles, like it is the case for 
nuclei or the electron gas, to mention just two prominent examples. 
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A Extension of linear response theory to unstable 



In the presence of an unstable mode of frequency lo + = iT the relations between response 
function have to be modified as compared to the common situation of a stable system. As 
we are going to show in this Appendix, some relations between response function and its 
dissipative part can be retained by making use of deformed paths of integration in the complex 
frequency plane, as encountered in the text for the contour C of (1161) . 

Causality requires that a physical quantity can be influenced by external fields only from 
their values in the past. The retarded function describing such a response must thus vanish 
for negative times. This may be achieved by an appropriate choice of the contour C + for that 
integral which defines the Fourier transformation from the frequency representation back to 
time: For negative times, it must be possible to close this integral in the complex frequency 
plane without enclosing any pole of the retarded response function x R (co>). This allows one to 
write: 



The contour C + could for instance be chosen parallel to the real axis if it only lies above the 
u + = iT (being the pole of x R having the largest imaginary part). In the case of the advanced 



modes 




(A.l) 
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response, in ( A.l ) all + signs are to be replaced by — signs, with the contour C_ lying below 
all poles of the advanced function. In the following we will write both versions in combined 
form, using an index X to distinguish between R(etarded) and the A(dvanced). Whenever 
the possible signs ± will appear, the upper one will relate to the retarded, the lower one to 
the advanced response. 

To obtain the equivalent of the well known relation between the frequency representations 



of response function and its dissipative part, one may start with the inversion of ( |A.l 



X X M = ±J dte iu4 2ie(±t)x"(t) \lmu\ > T , (A.2) 

which is restricted to frequencies whose imaginary part is larger (smaller) than T. In a second 
step, one may combine ( |A.2j ) for retarded and advanced response to find an expression for the 
dissipative part of the response function, 

*•<«> = s (L ^ e ~"* RM - L £ e ~*») ■ (A - 3) 

Under the assumption that retarded and advanced response do not share any common pole 
both contours C± may be deformed to a common one C without changing the value of each 
integral. One only has to ensure that C has all poles of the retarded response on one side and 
those of the advanced one on the other. For a stable system, C can be taken along the real 
axis and the formulas presented here turn to the ones known from the common literature. 
The final step is to insert the dissipative part ( |A.3| ) into ( |A.2j ) and to exchange the order 



of integration for time and frequency. Defining the analytic continuation (with respect to u) 
of the dissipative part of the response function by 

X» = ^(X»-X A H) (A.4) 

one obtains the familiar form 

x R(o,) = [ for i mu > o and uo ^ iT (A.5) 

JC 7T il — UJ 

This equation is the key relation for obtaining the integral representation ( CT ) for the fluctu- 
ation dissipation theorem from the analytic continuation of its representation (f45|). 



B The modified oscillator of the Drude regularization 



In ([46]) a frequency dependent friction coefficient was introduced to render some previously 
diverging integrals well behaved. The form chosen there, 7(0;) = 7(1— iu / vj d)^ 1 , is commonly 
associated to the Drude regularization. Replacing the friction coefficient in the oscillator 
response (|T2"D by this function the related dissipative part changes to 



//dm/A_ 1 W(7H+7*(W)) 



XoscV) 2M(cu 2 -u7 2 ) 2 + ^(7(cu)-7*H)(a; 2 -ro 2 )+7(cu)7*Hcu 2 (B ' 6) 
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It is apparent that with increasing uj the form on right hand side of ( |B.6|) drops faster than 
the one with constant friction. Therefore, the x"g C ru (u;) simulates better a correct behavior of 
the true response function at large uj, whereas the difference to the mere oscillator response 
function is small in the regions of maximal strength. 

The modified response function has one more pole as given by the new secular equation 



C - icu(— + 7 



Muj 2 + i uj 3 

ZU D 







Actually, its solutions depend on two parameters only, as can be seen by introducing x 
and rewriting ( (B.7Q to 



x 3 + x 2 
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This equation allows for an analytic solution, but it is more convenient to just solve it numer- 
ically. It is, however, instructive to evaluate them in first order in vj/wd, or, less lengthy, to 
solve directly ( |B.7[ ) perturbatively. Defining rj = r y/{2Mvj) and ( = vj/wd, one finds for the 
first two solutions 



uj 



(i) 



zu ( +>/signC — rj' 2 — irj + 2r/( 



^signC — rj 1 — i2r] 



uj 



(2) 
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It is easily recognized that for 77^ — > they turn into the solutions u± for the common 
oscillator, with uj± = ±JC/M — 7 2 /(2M 2 ) — ij/(2M). For the third solution one finds: 



UJ 



(3) 



-iw D (l - 2rj() 



(B.ll) 



A few remarks are in order here. First of all, it is interesting to note that it is not the ratio £ 
itself which matters but the combination with rj. It is intuitively clear that for larger damping 
the Drude frequency must be chosen larger. Secondly, we see that for a large wd, at least, 
the third pole lies on the negative imaginary axis. 

Let us finally just write down an expression for the full response function which serves as 
starting point for expressing the fluctuations in terms of the \E r — function: 



X 



dru 
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uj + imp 



M {uj - ujW){uj - ujW)(uj - uj®) 
where uj^ are the solutions of (p.Tp . 



(B.12) 
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